Quantum noise and entanglement generated by a local quantum quench 
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We examine the growth of entanglement under a quantum quench at point contacts of simple 
fractional quantum Hall fluids and its relation with the measurement of local observables. Recently 
Klich and Levitov proposed that the noise generated from a local quantum quench provides a 
measure of the entanglement entropy. Their methods were specific to non-interacting electrons and 
the generalization to interacting systems was left as an open question. In this work, we generalize 
their result to the Laughlin states. We investigate the noise generated in the current along the 
edge of a fractional quantum Hall state at filling factors v — 1/m, when a quantum point contact, 
initially closed, is fully opened at some initial time to = 0. We find that local quenching in these 
systems gives time dependent correlation functions that have universal behavior on sufficiently long 
time and length scales. We calculate the noise and full counting statistics for v = 1/m and find that 
in general, the entanglement entropy and noise generated are unrelated quantities. We also discuss 
a generalization of this problem to the critical quantum Ising spin chain. 
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I. INTRODUCTION 

Quantum impurity problems in 1+1 dimensions have 
long been the subject of intense study due to numerous 
applications in fields such as quantum wires, the quan- 
tum Hall effect and the Kondo problem. In particular, 
much research has been invested into understanding the 
behavior of particle transfer through an impurity and of 
growing interest experimentally and theoretically are the 
statistics of the fluctuations in the transferred charge. 
Most efforts have been concentrated on the shot noise 
which is sensitive to the quantization of charge and can 
give useful information about the systemii*2i2i^i5i& In re- 
cent years, it was realized that the full counting statistics, 
probability distribution P(q) of the transmitted charge q 
in a given time window At, contains more information. 7 
Along with the shot noise, the full counting statistics also 
contains information about the higher order correlations 
which are of interest experimentally. It has been sug- 
gested that the third moment may give a more reliable 
measure of the charge than the shot noised Recent ex- 
perimental efforts suggest that the higher moments of the 
full counting statistics may be accessible^ ' 10 ' 11 ' 12 

In cases where Fermi liquid theory applies, the full 
counting statistics has been studied extensively (see Ref. 
[l3| for a review). Systems of strongly interacting elec- 
trons, though of importance in 1+1 dimensional quan- 
tum wires, the quantum Hall effect and Kondo problem 
have been less well studied. Using advanced methods 
such as the thermodynamic Bethe ansatz, for the Laugh- 
lin states exact results for the charge current through 
an impurity^ as well as the noise^ the heat current?^ 
and even attempts at the full counting statistics were 
reportedji 7 - Exploiting the power of boundary conformal 
field theory, the full counting statistics in double quan- 
tum dots was studied^ 

A subject of recent interest is the behavior of quan- 
tum entanglement, and in particular of the entanglement 



entropy, in condensed matter systems (for a recent re- 
view see Ref. [Til ■) The entanglement entropy has been 
shown to exhibit universal scaling behavior near quantum 
critical points. 20 ' 21 ' 22 ! 23 ' 24 Universal scaling behavior of 
the entanglement entropy is also present in topological 
phases of matter, such as fractional quantum Hall fluids 
and their generalizations, as well as in the related topo- 
logical quantum field theories. 25 ' 26 ' 27 However, the entan- 
glement entropy of a macroscopic quantum system is a 
highly non-local quantity which is difficult to measure. It 
has remained a challenge to find an experimentally viable 
protocol to measure the entanglement entropy. 

Several recent results have suggested that the behav- 
ior of point contacts in quantum critical systems, and 
topological phases, may offer a way to measure the en- 
tanglement entropy. Fendley et al have shown 2 ^ that the 
change in the entanglement entropy of topological FQH 
fluids at a point of constriction is related (in fact the 
same) to the change of the Affleck-Ludwig entropy^ 9 - of 
the coupled edge states of the FQH fluid at the point 
contact. However, the Affleck-Ludwig entropy itself is 
difficult to measure. 

More recently, Klich and Levitov^ have shown that, 
at least for a system of free fermions, it may be possible 
to measure the growth (in time) of the entanglement en- 
tropy upon a quantum quench by monitoring the noise in 
the charge current through the contact. Specifically, they 
proposed that the second cumulant of the full counting 
statistics is related to the entanglement entropy. They 
considered the following protocol: first the two subsys- 
tems described by non-interacting fermions lying on ei- 
ther side of an infinite strength impurity are completely 
decoupled. At time to = 0, the impurity is removed and 
the subsystems are allowed to exchange particles. Finally 
at time t\ , the impurity is reinstated and blocks the flow 
of particles. They find that the shot noise generated is 



2 



given by 

1 At 

Snoise =o lo S — I 1 - 1 ) 
o r 

where r is some short time cutoff. 

Such a protocol is a suggestive procedure for measur- 
ing the quantum entanglement between two halves of a 
system. In the initial state the two subsystems lying on 
either side of the impurity are completely independent, 
but following the quench, quantum entanglement is dy- 
namically built up by the exchange of particles with van- 
ishing net flow. A measure of the entanglement between 
the two subsystems is provided by the entanglement en- 
tropy. The entanglement entropy for such a scenario was 
recently calculated using conformal field theory (CFT) 
and shown to increase logarithmically with time, with a 
universal coefficient proportional to the central charge c 
characterizing the CFT^^^ 

c At 

S ent = |log— (1.2) 
3 a 

where a is a short distance cutoff. By comparing Eqns. 
II. HI. 21 and noting that for non-interacting fermions, the 
central charge is c = 1, it is quite suggestive that the 
noise is a measurement of the entanglement. However, 
it should be noted that while the methods of Calabrese 
and Cardy are general, the results of Klich and Levi- 
tov are specific to non-interacting electrons. The degree 
of general validity (e.g. for interacting systems) of the 
Klich-Levitov protocol is presently an open question. 

Recently there has also been a surge of interest in 
quantum quenching, i.e. where the parameters that de- 
scribe the dynamics of the system are changed over a 
short period of time either locally or globally. In cold 
atoms, it has been seen experimentally that quantum 
quenches may show interesting features i 34 ! 35 ' 36 Theoreti- 
cally, quantum quenching has also attracted a good deal 
of attention in recent yearsi 31 i 37 ' 38 i 39 ' 40 It was found re- 
cently that in a global quench, where the eigenstate state 
\tpo) of a Hamiltonian Hq is evolved by a different Hamil- 
tonian H, the correlation functions were found to display 
universal behavior characteristic of the Hamiltonian H if 
it was tuned near criticality4i 

In this work, we study a local quantum quench where 
the point contact between two v — 1 jm fractional quan- 
tum Hall states is instantaneously opened (see Fig. []}. 
We examine the noise and the full counting statistics and 
extend the results of Klich and Levitov[3(| to a system of 
interacting electrons with a dynamic impurity. No exter- 
nal bias is applied. Here the bare strength of the impurity 
coupling undergoes a sudden change between two values, 
the first corresponding to fully reflecting boundary con- 
ditions and the second to fully transmitting boundary 
conditions. The quantum point contact (QPC) generates 
an effective impurity in the Luttinger liquid description 
of the edge, thus allowing complete control over the im- 
purity strength by tuning the voltage applied to a side 



gate. Such a system is therefore a promising candidate 
for an experimental exploration of quantum quenches. 

Our main result is that the second moment of the full 
counting statistics, and therefore the noise, has a simi- 
lar form to the entanglement entropy. This result holds 
even for interacting (Luttinger) systems. However, we 
also find that this correspondence appears to be coin- 
cidental since, in general, in addition to the noted de- 
pendence of the central charge, the noise depends also 
on other universal quantities of the underlying confor- 
mal field theory of the system. There is however a con- 
ceptual difference between static and dynamic entangle- 
ment. Static entanglement refers to the measure of the 
non-local correlations associated with the observation of 
a part of the system without disturbing it and is en- 
coded in the von Neumann entanglement entropy. On 
the other hand, when a physical system is changed, as in 
a quantum quench, the ensuing dynamical entanglement 
measures the non-local correlations that develop upon its 
time evolution in the quenched system. (See, e.g. Ref. 
[52]). At the computational level, the main method un- 
derlying our results is the use of a boundary condition 
changing operator to generate the transition between the 
fully closed QPC and the fully open one. This transition 
is related to a change in the sign of the odd-boson density 
operator going through the QPC, thus is marked by the 
appearance of an orbifold theory 

The method in its most general form will be described 
in section IIVI In sections IIIII and IIVI we describe the 
method for v = 1/2 and 1/3. For the particular case 
of v = 1/2, the result can be derived explicitly using a 
Majorana fermion representation for the re-fermionized 
Hamiltonian. Finally, in section IIVBI we calculate the 
noise, and in IIV CI the full counting statistics. In Section 
IVl we sketch the calculation of the noise of the energy- 
momentum in the quantum Ising chain, showing that in 
contrast to the entanglement entropy in a 1+1 dimen- 
sional critical system, the noise does not always grow 
logarithmically. The details of the calculations for the 
Laughlin v = 1/2 (bosonic) state are presented in Ap- 
pendix [X] The modular S'-matrix and fusion rules for 
the orbifold CFT (used in Section [TV]) are given in Ap- 
pendix [Bj The relation between the Schwinger-Keldysh 
formalism and entanglement entropy is summarized in 
Appendix [CJ 

II. DESCRIPTION OF THE MODEL 

The Lagrangian describing the edge of the quantum 
Hall effect for the Laughlin states at filling factors v — 
1 jm can be written in the terms of left and right moving 
bosonic fields, (f> L / R , in the form 

c = ^-d x <p L (d t - d x )cb L - ^-d x <p R (d t + d x )4> R (2.1) 

where x £ [— £,£], and the velocity v of the edge was 
set to one. For decoupled edges, the physical Hilbert 
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FIG. 1: Single quench at t = 0. (I) Fully reflecting for t < 
(A — > oo), (II) Fully transmitting for t > (A -» 0). 



space of each chiral boson is generated by j, the £7(1) cur- 
rent and the vertex operator e % ^^ L i R (see Ref . |43ll44| ) . 
This vertex operator generates the charged excited states 
of the theory which are fundamentally generated by 
adding or subtracting an electron from either the left 
or right edge i 43 ' 45 Requiring that the charged excita- 
tions generated by e l ^^ L ' R are bosons (electrons) im- 
plies that they satisfy the (anti-) commutation relation 
sets g = to an even (odd) integer. The requirement that 
the charge of the operator is precisely e forces the rela- 
tion v = l/m. Representations of the theory are local 
with respect to the U(l) current and charged excitation 



L/R 



e *^/g<pL/R w hich is identified with the electron 



(g = m odd) or boson (g = m even). 

The QPC introduces backscattering at a; = 0, with 
the most relevant term being the quasi-particlc tunneling 
term 



C! = X6(x) cos [V"((t> R - L )] 



(2-2) 



where A is the tunneling strength. The charge densities 
at the left and right edges are respectively 



PL 



PR 



0, 



(2.3) 



Note that <\>i,{x + 21) = 4>l{x) + 2tt^/vQ l , so in partic- 
ular the electron operator is well defined provided Q L is 
an integer; hence Ql carries the meaning of the number 
of quasi-particles on the edge. In units of the electron 
charge e, the charge carried by the edge is vQl- 

We now perform the following transformation into the 
odd and even basis 



1 



<t>o(x,t) = — = (<j>L(x,t) ~ (f>R(-X,t)) 



<j>e(x,t) = -j= (<j>L(x,t) +(j) R {-X,t)) 



(2.4) 



the even boson <p e decouples from the impurity, and we 
ignore it in the rest of the paper. Note that 4> / e are 
both left-moving. The odd and even charges are related 



to the original L/i?-charges through \f2Q° = (Ql — Qr), 

V2Q e = (Ql+Qr)- 

The experimental setting we wish to consider is one 
where the two edges of the quantum Hall liquid are sep- 
arated initially. No bias is applied and the temperature 
is taken to be zero. At some time, the two halves of 
the system are (suddenly) connected and we wish to con- 
sider the statistics of the charge transferred from one side 
to the other upon this change, i.e. a quantum quench. 
The system is then evolved with the Hamiltonian with 
a Hilbert space described by states with transmitting 
boundary conditions at the quantum point contact. 

Let P(q) be the probability that charge q is transmit- 
ted through the point contact in time At. Then one 
defines the generating function x(A) = X^-oo P{l) e%Xq 
which encodes all moments (or cumulants) C m of the 
probability distribution P(q), 



C m = (-idx) m log X (X) 



(2.5) 



A=0 



In particular, the second moment Ci is related to the cur- 
rent fluctuations. The generating function can be written 
as 



xW = ({ 



e i\q(At) e -i\q(0) 



}) 



(2.6) 



where the operators are ordered on the Schwinger- 
Keldysh contour?! This can also be written as a trace 
over a complete set of states. 



X(A) = Tr {p U^ Xq Ue- iXq ) 



(2.7) 



where po is the initial density matrix, and U is the time 
evolution operator for time At after the point contact 
is opened. In Eq. (|2.7p it should be understood that the 
trace is taken over the initial states. 

If the scattering time at the quantum point contact is 
short compared with the entire time evolution, then it 
was shown that xW can be written asi^ 

X(A) = det (1 + n{S^e lXq Se- lXq - 1)) (2.8) 

where n is the initial distribution of states and S is the 
reflection-transmission matrix at the quantum point con- 
tact. In Ref.[3(|, it was found that for non-interacting 
electrons, the generating function is given by 



X(A) 



= e ~ — C 2 



(2.9) 



where C 2 = £ log (f). 

In this work, we generalize this protocol to the v = 
l/m Laughlin states. We note that while the full count- 
ing statistics has been found in many other strongly cor- 
related systems with a static impurity via the thermo- 
dynamic Bethe ansatz ) 14 ' 15 i 17 ' 47 the application of this 
method to finding time and space dependent correlation 
functions is in general difficult and remains largely not 
understood. 
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We take a different approach in tackling this problem. 
This relies on the realization that one should consider an 
extended Hilbert space by introducing a boundary chang- 
ing operator to the edge theory and that physical charged 
excitations in the initial Hilbert space are those conserv- 
ing total charge. In the next section, we examine the 
noise in v — 1/2 which can be refermionized and solved 
explicitly. We note an interesting structure which we 
then extend to the other Laughlin states at v = 1/m. In 
doing so, we calculate the noise and full counting statis- 
tics for the Laughlin states. 



III. NOISE AT FILLING FACTOR 1/2 

We start our discussion with an exactly solvable case, 
the bosonic quantum Hall effect at filling factor v = 1/2. 
Our strategy is as follows: re-fcrmionizing, we first re- 
cast the problem in terms of a quadratic fermionic ac- 
tion. Then, by writing the fermionic field in terms of two 
real Majorana fermion fields, we find that only one of the 
Majorana fermions interacts with the impurity. The two 
fixed points associated with the fully transmitting and 
the fully reflecting impurity, translate into anti-periodic 
and periodic boundary conditions imposed on this Ma- 
jorana fermion. This leads us to identify the operator 
that takes the system between the two fixed points as 
the a operator of the chiral Ising model. Armed with 
this knowledge, we then proceed to use this formalism 
to calculate the noise for the case of a quantum point 
contact which is suddenly opened at time t = 0. 



A. Majorana fermion representation 

The odd boson theory can be re-fermionized in terms of 
a fermion field ip (see Appendix [X| leading to the Hamil- 
tonian 

H a = J dx ip^ (x)id x ip(x) 

+V2nXS(x) [^(x) 7 + 7?/' t (x)] (3.1) 

The crucial step is that the Hamiltonian can be written 
in a Majorana representation, ip = (r/i — ir\?)l2 where 
•qi and 772 satisfy {^(x), 77,(2')} = 2SijS(x — x'). In this 
representation, only 772 interacts with the impurity, as is 
clear by inspecting the equations of motion, 



Integrating the equations leads to the following boundary 
condition across the impurity, lob = 47rA 2 , 



id t ir]\ = id x ir]i 

idtir]2 = id x ir\i + 2V2ir\j5(x) 
ida = 2%/2^AiT? 2 (0) 



%(0 + ) 



LJ B 



(3.3) 



11a — lob 

This boundary condition is continuously interpolating 
between the fully reflecting and fully transmitting bound- 
ary conditions for the odd boson 



LO B —> OO ?72(0 + ) = -772(0 ) 
LOB^O 772(0+) = 772(0-) 



(3.4) 



(3.2) 



Similarly, at x = i, we have the boundary condition, 



= (3.5) 



The problem is thus mapped to a single chiral Majo- 
rana fermion on a circle of circumference 2£. For the 
reflecting case we have anti-periodic (AP) boundary con- 
ditions, while for the transmitting case, we have periodic 
(P) boundary conditions. 

The free Majorana field theory contains altogether 
three primary fields: the identity field 1, the Majorana 
fermion field 77, and a a operator (known as the Ising twist 
field) creating a branch cut on the Majorana fermion 77. 
As before, we shall label the two copies of fields present 
here by an index i = 1, 2: 1,-, r]i, and Cj. To go from pe- 
riodic to anti-periodic boundary conditions on 772 , we use 
the (T2 operator to create a branch cut on the fermion: 
we imagine that the edge Majorana fermion encloses a 
"bulk" in which 02 operators can be introduced, appro- 
priately changing the boundary conditions on the Majo- 
rana fermion. We thus identify the operator taking the 
system between the two fixed points associated with the 
QPC as the a% operator in the 772 Majorana fermion field 
theory. 

In particular, to go from anti-periodic to periodic 
boundary condition, introduce two 02 operators in the 
"bulk" surrounded by the 772 edge (see Fig. [5]). One of 
them then approaches the edge at x = and a bulk- 
edge coupling is introduced. This coupling is mediated 
by virtual processes involving the tunneling of topolog- 
ically trivial modes between the bulk and the edge^. 
Finally, when A — > 00, the edge circumvents the operator 
completely, and switches to AP boundary conditions. In 
the sudden approximation, we can imagine two spatially- 
separated CT2 operators introduced simultaneously: one 
inside the edge and the other outside. 



B. Calculation of the noise the noise produced in that situation. 

We assume that the QPC is initially closed, then at 
time t = it is fully opened, and proceed to calculate 



FIG. 2: Boundary changing operator - fully transmitting to fully reflecting. Two 02 operators are drawn from the vacuum. 
Tunneling (of strength A) is introduced between one of the 02 operators and the edge. Finally, in the limit A — > 00, the edge 
circumvents the a% operator. 



The current operator through the junction can be de- 
fined as the rate of transfer of charge between the two 
edges, measured by the odd boson 



I = d t po = i[H, po] 



(3.6) 



where p = ift'ip. Using the equations of motion, Eqns. 
(|3.2p . we can write the current in terms of the Majorana 
fcrmions as 



In the limits, ub —* 00 and — > 0, we get 



% (0)772(0+) 



t < 
t > 



(3.7) 



(3.8) 



where t = time of the quench, and I b is what we call 
the backscattered current. 

In the following, we shall define the charge transmitted 
through the junction during a time window starting at 
t = —00 and ending at some arbitrary positive time At, 



Al 



dtl b {t) 



Al 



dtl b (t) 



(3.9) 



Since there is no bias, the net transferred charge is zero, 
(<?At) = 0. However, there will be noise due to the open- 
ing of the QPC. The noise is given by 



(Qlt) 



Al 



dt-i 



Al 



dh Mowo)> (3 ' 10) 



Besides its role in establishing a non-zero current opera- 
tor, the operator 02 also changes the boundary conditions 
on 772. However, these boundary conditions {i.e. the ef- 
fects of the operator cr 2 ) will affect 772 only after a long 
time-scale set by £/v. In the limit I — > 00 we can safely 
neglect the effects of these boundary condition changing 
operators. Upon the introduction of a UV cutoff S, the 



short-time switching scale needed to regulate the inte- 
gral, the second cumulant C 2 can be obtained. 



(Qlt 



2 ~ 2 I *i / ^y— —-—2 (3.H) 



2^ 2 Jo Jo (S + i(ti-t 2 )) 
The expression can be easily integrated to get 



(Q 2 At) = 7TT 1o ' 



v , At 2 + 5 2 



In the limit of large At this reduces to 



Co 



v , At 
— lo g-T 

7T^ 



(3.12) 



(3.13) 



The essential point is that the theory of the two edge 
states coupled by the local QPC must have a local bound- 
ary changing operator 02 that takes the system from one 
fixed point to the other. In particular, the operator <r 2 
maps the ground state of the system with a closed QPC 
to an energy eigenstate of the same system but with the 
new boundary conditions, an open QPC. The ensuing 
growth of the entanglement entropy is due to the evolu- 
tion of correlations of the excitations in this state. 42 



C. Bosonic Picture 

The previous discussion can also be understood in 
terms of the boson. This picture lends itself to generaliza- 
tion to the other Laughlin states which are also described 
by free chiral bosonic theories. Here, we show that from 
the bosonic picture, one can deduce that a operator must 
be introduced for one of the Majorana fermions to change 
the boundary condition from anti-periodic to periodic. 
We then argue that a similar boundary condition chang- 
ing operator can be included in the bosonic theory. 

One can understand the relationship between boson 
and fermion boundary conditions and the nature of the 
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twist operator by examining the partition functions. Spe- 
cializing to v = 1/2, the Lagrangian becomes 

C = -^d x Mdt ~ d x )cf> + \8{xy^ + h.c. (3.14) 

4-7T 

The boundary condition at the impurity site changes 
from A — ► oo perfectly reflecting (Neumann boundary 
conditions) to A — > perfectly transmitting (Dirichlet 
boundary conditions). The Neumann boundary condi- 
tion simply fixes the field <j) to be continuous across the 
boundary while the Dirichlet boundary condition allows 
the field to be discontinuous. 

Imposing Neumann boundary conditions at the ends 
of the system for all times, the mode expansion for the 
odd boson can be found. For Neumann (N) and Dirichlet 
(D) boundary conditions at the origin respectively, 



(a) 



(po(x,t) = 0o + —Q°(t-x) 



n£Z 



t> (x,t)=iJ- 

V rGZ+i 



i(t— x) 



(TV) 



(3.15) 



The partition function in each sector can then be com- 
puted for each set of boundary conditions and one finds 



A 



7 _ 03(g) _ 7 f 7 f 



(3.16) 



A^O 



7 UlWsjq) _ f f 

/jNU ~ K 1 ^27Z\ - ^PA^AA 



r] 2 (q) 



where Z- are partition functions for a single c = i holo- 
morphic fermion with boundary conditions i in time and 
j in space, A anti-periodic and P periodic. To go from 
A — > oo to A — > 0, we see that exactly one Majorana 
fermion changes anti-periodic to periodic boundary con- 
ditions in time. Now, the partition function Z pA is also 
related to the partition function Z AA by an insertion of 
the fermion counting operator (—1) in the time direc- 
tion. By a conformal transformation, this operator is 
related to the spin operator of the Ising model inserted 
in the space direction for one of the Majorana fermions 
that compose the theory4^ This confirms the results of 
the previous section. 

While this identifies the boundary changing operator 
as the spin operator in the Ising model, this is a spe- 
cial feature of v = 1/2 state. At other filling fractions 
the boson can no longer be fermionized, but what is 
independent of the filling is the boson mode expansion 
and the partition functions Znd, Z^n- To generalize 
the calculation of the noise beyond v = 1/2, one should 
identify the boundary condition changing operator in the 



t=0 



(b) 



FIG. 3: (a) The Schwinger-Keldysh contour can be regarded 
as the time contour of a conformal field theory by an analytic 
continuation, (b) The path integral can be thought of as 
over the entire surface parameterized by (a;, r) with changing 
boundary conditions at t\ and 75. Details in Appendix [Cl and 

nam. 



bosonic theory that takes Neumann boundary conditions 
to Dirichlet boundary conditions. 

A priori, it may seem strange that a c = 1 free bosonic 
theory must include a boundary changing operator, but 
the necessity of this operator can be seen at the level of 
the path integral. The path integral should be time or- 
dered on the Schwinger-Keldysh contour, but by an ana- 
lytic continuation, one can regard the Schwinger-Keldysh 
contour as the time surface for a conformal field theory 40 
(see Figj3|). In going from the density matrix to the 
path integral formulation, one must consider an extended 
Hilbert space of the Hamiltonian TL which is spanned by 
states with all possible boundary conditions. The bound- 
ary changing operators appear as unitary operators that 
map whole sectors of the extended Hilbert space with 
particular boundary conditions to sectors with different 
boundary conditions. 50 Details of this analytic continua- 
tion can also be found in Appendix [Cl and in Refs. 40lf4ll |. 

Folding the system on the (x, r)-world sheet, the local 
quenching problem can be regarded as a boundary con- 
formal field theory with a changing boundary condition, 
Neumann to Dirichlet at the origin. Such a boundary 
conformal field theory was studied in the context of the 
X-ray edge singularity and quantum wiresi 50 ! 51 The Neu- 
mann boundary condition implies that the odd boson is 
continuous at x = while the Dirichlet boundary condi- 
tion implies that the odd boson is discontinuous at x = 0. 
As discussed by Affleck and Ludwig^ this can be written 
as: 



o (O + ,t<O) = o (O-,i <0) 



o (O+,i>O) = -0 o (O-,i>O) 



(3.17) 



This has an interpretation of a branch cut running on 
the (x, r)-world sheet across which the odd boson is dis- 
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continuous. In Ref. 50] it was noted that this branch 
cut corresponds to an insertion of a dimension h = 
operator at the quench times, t\,t%. 

Observing that the boundary conditions are indepen- 
dent of the radius and following arguments by Ref. [50( , it 
is suggestive that there is dimension h = twist oper- 
ator in the odd boson theory (independent of the boson 
radius) which changes the boundary condition from Neu- 
mann to Dirichlet. Given such generality, one wonders if 
there is a larger algebraic structure describing the local 
quench which includes the boundary condition changing 
operator in its operator algebra. In the sequel, we show 
that this structure is described by an orbifold conformal 
field theory. 



can be found in any standard text 







-19.52 



h = —, 

4m' 



where 



= e ±n 



2m0 and 



(4.1) 



Orbifolds 



IV. EXTENSION TO OTHER LAUGHLIN 
STATES 



Much of the previous discussion was limited to the v = 
1/2 case and was special in that it could be formulated in 
terms of Majorana fermions. We argued previously that 
to generalize the previous results to the other Laughlin 
states, one needs to understand the theory in terms of free 
chiral bosons and to work with an extended Hilbert space 
where a dimension h = boundary condition changing 
operator is included. Such a Hilbert space is described by 
an orbifold conformal field theory, but before delving into 
orbifold theories, we need to describe the initial Hilbert 
space of the system which we will extend in the next 
section. 

The total charge of the system is conserved in the 
quenching process so that physical charged excited states 
are those preserving the total charge. Initially, one 
knows that in the limit of two decoupled edge states, 
the charged excitations that conserve total charge are 
processes where an electron tunnels from one edge to the 



other, d 



-,T,t 



^ r ^B + h.c. In the even and odd basis, one 



can identify the vertex operator e ± l V2g<P° with g = m as 
generating these excited states. By charge conservation, 
the even sector does not have charged excited states, and 
its Hilbert space is generated by the U(l) current. How- 
ever, charged excited states are allowed in the odd sector 
by electron tunneling, so the physical Hilbert space is 
generated by the U(l) current as well as the vertex op- 
erator e ±i%/3 ™*° . 

Representations of this Hilbert space can be con- 
structed with the requirement that they be local with 



respect to the generators j 



and V± = e 



This construction is well known and extends the U(l) 
Kac-Moody (KM) algebra to conformal field theories de- 
scribed by the A m -series. The primaries in the theory 



The twist operator a exists explicitly in the fermion 
theory as spin operators, but by considering a Z2 orb- 
ifold of the free Gaussian theory, a twist operator can be 
included in the boson theory. On the (x, r)-world sheet of 
the odd boson, one has that the <p (z) = —<p (z) as z ap- 
proaches the branch cut running between the two quench 
times. Imposing 4> o (0 + ) = — </> o (0~) on the (x,r)-sheet 
has consequences for the target spaced For a compact 
boson, <fio E S 1 = [0,2irR], this means that the target 
space is identified under the action of the discrete group 
%2- The theory must include an operator which takes 
4>o — > ~<t>o- 

As argued before, the theory describing the odd boson 
is a conformal field theory in the A m -series. To this the- 
ory, we wish to allow for an anti-periodic boson and need 
to identify the target space under the group Z2. One 
needs to consider the orbifold j4 m /Z2. The partition 
function for orbifold theories under the group Z^v was 
found by Ginsparg. 54 Decomposing this partition func- 
tion into A m characters and observing the transforma- 
tion properties under the modular S : r — > — 1/r and 
modular T : t — > r + 1 transformation allows one to de- 
duce the primaries, the modular iS-matrix of the theory 
and fusion rules. 52 ' 55 The primaries of the theory are, 

1 <-> h = Q 



1,2 



h= — , fc=l, ..JV-1 
4m 



01,2 



n,2 



16 
9 

16 



(4.2) 



The Z2 transformation g acts on these representations 
as g : [cf>k] — > [4>2m—k] an d so there are two fixed points 



8 



[<j>o] and [4> m ]- The operators <fim correspond to a sector 
invariant under (j) — * ~4> (i- e - 't'm — cos(y / 2m<^ ) )and 
a sector that is broken by the transformation <p Q — > — O 

a( 2 ) 



(i.e. 0m = sin(v 2m0 o ) ). Physically, one can under- 
stand them as the tunneling and current operators for 
the odd boson respectively. The operator 4>k are in- 
variant under the transformation <ftk — cos( , fc 



55 



cti.2 correspond to the trivial and non-trivial represen- 
tations of Z2. For m even, it turns out that the <Ji 
fields are self conjugate while for m odd, o\ and 02 
are each others conjugates^ Using the fusion rules (see 
Appendix [B]) and decomposing the partition function 
Zbd in terms of the characters of the A m j%2 theory, 
it can be seen that the operators act as boundary 
changing operators. For m even, fusing with o\ gives 
Zdd — Znd —f ai Znn while for m odd one has 

ZdD — ►o-i Znd -^<t 2 Zmn^ 

For v — 1/2, the orbifold theory is A2/Z2 which is 
isomorphic to the tensor product of two c = i Vira- 
soro algebras (i.e. two decoupled Ising models). 55 This 
was exactly the behavior observed before. The partition 
function for the odd boson one observed could be decom- 
posed into a product of c = i partition functions. In 
addition to this, the primaries of the theory are 

1 <-> h = 



in the non-trivial representation and non-local with re- 
spect to the fermion 1J2 (see Fig. [5]). When the interac- 
tion strength is switched, precisely one spin field inter- 
acts with either 7712 giving a twisted boundary condition. 
This behavior was also seen in the odd boson partition 
function and explicit calculation in Section [IIII 

For v = 1/3, the theory is described by the orbifold 
theory A3/Z2 which corresponds to the Z4 parafermion 
theory^ This can also be seen by comparing the fusion 
algebra and operator content of the two theories: 



ft = 



01.2 



= 1 

_ 3 
~ 4 

1 

~ 12 
1 

~ 3 
1 

" 16 



0"1,2 



n,2 



h = 1 

h= 1 - 
2 



ft=l 
16 

16 



(4.3) 



We see that the algebra naturally includes the twist fields 
0X2 and t\ 2 and that there is an operator corresponding 
to the quasiparticle tunneling operator, 



3 X 3 = 1 



= 1 



,< (!) 



4 2) 



3 



o\ x a 2 



4« 



(4.4) 



The operator <j)x has a fusion rule consistent with its in- 
terpretation as cos(|0 o )i^ 

In addition, we see that the dimension h = g field we 
can regard as the dimension of the two spin field com- 
posite. Hence, one of the er's originally at the center can 
be taken to be in the trivial representation <j\ of Z2 non- 
local with respect to the fermion field r\\ and the other 



T l,2 



h= — 



9 

16 



(4.5) 



The fusion rules here arc somewhat different than for 



even. The operator algebra of 1, j and yj 3 
symmetry, and it is given by 



now has a Z4 



3 X 3 = 1 



f x $ 



3 4 X) x <4 2) = 1 



ai X (72 



1 -I 



(4.6) 



The fusion algebra for the <j>k vertex operators is 
unchanged^ 5 - The fusion rules and modular S'-matrices 
for other filling fractions are given in general in Appendix 

m 



B. Calculation of the Noise 

One can then go about computing the noise in the 
backscattering current in the language of bosons. In the 
original left-right basis, the backscattered current is given 
by 



j t An 1 J,n\ 1 out , „out \ 

*b - (Pl +Pr)~ (PL + PR ) 



(4.7) 



pout are chgugg densities measured at a position after the 
point contact while p] n are charged densities measured 
at a position before the point contact. In the rotated 
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even-odd basis, the backscattering current can be written 

as, lb — : ^-{d x 4>°_ — 9x0+) where <p± denotes after and 
before the point contact at x = 0. To connect with the 
boundary conformal field theory formalism, one can fold 
the system and regard dcj)(x < 0, t) as a right moving 
wave B(f>(x > 0,t). The backscattering current is then 



J 



h = \Zj^(j° ~ 3°) where x > now. This form has the 
advantage of explicitly involving a primary field of the 
theory. Intuitively, the backscattered charge is simply 
the difference between the (incoming) left moving charge 
density and the (outgoing) right moving charge density. 
The noise is, 



({h(r 2 ),I b (r 1 )}) a 



(.f(lT 1+ 6 e )3°(lT 2 

(j (iT 1+ S e )]°^r 2 - 



-*«)>« 
I — 



(?( l r 1 -5 e )j°(^2 + ( 5 e )) a + (r 1 ^r 2 ) 



(4.8) 



where (•••)„ denotes a correlation function with boundary 
condition a at x — (Neumann/fixed or Dirichlet/free) 
and fixed (Neumann) at x = i. Here n = t\ — ie and 
t 2 = t 2 + it when t\ is on the top contour and t 2 on the 
bottom contour and vice versa in the second contribution. 
S e is a short distance cutoff, 6 e — » 0. In the presence of 
a boundary condition a these correlation functions have 
been computed using methods of conformal invariancc 
and their relationship to the modular S'-matrix of the 
bulk theory were found in Refs.[57l)58j|. 



{<j>{z l )<j>{z 2 )) a = M*iM*2)> 




for xy > 
for xy < 

(4.9) 

where z\ = ir±+x and z 2 = ir 2 +y and <p is a primary field 
with scaling dimension A. By knowing which primary 
field corresponds to boundary condition a one can use 
Cardy's formula to compute the correlation function in 
the presence of the boundary. 

For v — 1/2, we showed that the theory can be thought 
of as A2/Z2 which is isomorphic to a theory of two de- 
coupled Ising models. The boundary state correspond- 
ing to free boundary conditions (Dirichlet) is associated 

I 



with the primary a while the fixed boundary conditions 
(Neumann) is associated with the identity or e primary 
operators. Using the 5*-matrix for the Ising model, its 
easy to see that the correlation functions with boundary 
condition a have the prefactors, 



fixed/N: 
free/D: 



S° /S° ~ 



(4.10) 



For other filling fractions, the identification of boundary 
states and boundary conditions is more complicated but 
we find that a similar structure is present. 

Since the characters for the A m /1 2 theory are known, 
the partition functions Zjyjy and Z^p can be decom- 
posed into a sum of characters from which one may iden- 
tify the boundary states i 56 ' 59 Note that here we worked 
with the T-dual (r -> §) of the theory in Refs. [Hil] 
so that D <-> N. For m even, Neumann and Dirich- 
let boundary conditions, it was found that the orbifold 
boundary states are, 



I AT) 



|0) + li) + 7i( 



\D) 



m|0) 



m / 
~2 V 



rn — 1 



k=l 



(4.11) 



For m odd, the analogous states can be written down 
using the appropriate S- matrix (see Appendix |B|) . From 
the modular S'-matrix, one then can identify the Neu- 
mann boundary condition with the highest weight state 
</>fc and the Dirichlet (free) boundary condition with the 
a field. Note that the boundary states above are not 



the usual boundary states but are Z2 invariant linear 
combinations^ However, the fact remains that they 
transform in the [4>k] and [a] representations respectively. 
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For v — 1 jm one finds that it is generally true that 



fixed/N: 
free/D: 



For arbitrary v, one finds that in the limit e — > 0,S e 
the noise from the local quench is given by 



.12) 



0, 



{{I b (r 2 ),I b (ri)}) 



l 



2tt 2 (S + i(h - t 2 )) 2 



(4.13) 



in agreement with the calculation for v = 1/2 so that the 
second cumulant has a similar form as Q3.1ip . 

A non-trivial check of these results can also be done by 
considering the boundary condition on the odd currents. 
In terms of the odd currents, the reflecting and trans- 
mitting boundary condition were given by j° — j° and 
j° = —j° respectively^ Imposing either of these bound- 
ary conditions, one sees that the left-left and right-right 
correlation functions are unaffected but the left-right cor- 
relation function changes by a minus sign. 



l l 



2 (*i - z 2 y 



(4.14) 



2 (zi - z 2 ) 



so that when the point contact is closed, the noise mea- 
sured is zero while when it is opened, the noise is non- 
zero. This is in agreement with the results using the 
operator algebra of the A m /Z 2 orbifold. Note that this 
correlation function is special in that the boundary con- 
dition involves the primary fields j° and j° that appear 
in the desired correlation function but Eqn. (|4.9|) is more 
general than this. 

While the second cumulant can easily be shown to grow 
logarithmically in time, our result shows that the coeffi- 
cient of this logarithm depends on the modular S-matrix 
elements of the theory and not just to its central charge, 
so that even for c = 1 the second cumulant does not fully 
give the entanglement entropy. 



even in the interacting case, there exists a simple relation 
between the FCS and the entanglement entropy. 

To avoid boundary effects, it is convenient for our pur- 
poses to consider a Corbino (annular) geometry, with a 
single QPC inducing backscattering between the edges. 
In this new geometry, we should actually reverse the 
protocol appearing in previous sections: first the QPC 
is open, then is quickly closed, and remains in this 
state during a time period At, at the end of which it 
is opened again. Going to the odd/even basis for the 
two edge bosons, the charge carried by the odd boson, 
Qo = f dxp = (4> o (0 + ) - O (O~)), is controlled by 
the singularity induced by the branch cut for the closed 
QPC (assumed to be located at x = 0). The full count- 
ing statistics is then given by the following correlation 
function, 



xW = {{ 



e iXQ„(At) ^ e -i\Q o (0) 



}) 



(4.15) 



which, using the observation above, can be reduced to 
the calculation of vertex operators at the vicinity of the 
QPC. This leads to the main result, 



X(A) = exp 



2 2tt 2 



log 



At 2 



5 2 



(4.16) 



which is the typical generating function for a Gaussian 
distribution. Consequently, only the second cumulant 
(i.e. the noise) is non-zero. Interestingly, except for the 
explicit dependence on the filling factor v, the full count- 
ing statistics has the same form as the non-interacting 
electron case. 

Writing the boson in the even and odd basis and fold- 
ing the system, this problem can be mapped onto the 
boundary sine-Gordon model. In the basis of kink and 
anti-kink solitons of charge q and —q respectively, dis- 
tribution functions can be found via the thermodynamic 
Bethe ansatz for the entire range of Aj 15 ' 17 The prob- 
lem is then similar to the non-interacting electron prob- 
lem only the particles are collective modes rather than 
the original electrons. However, in this picture, the full 
counting statistics is difficult to compute. The reflection 
and transmission coefficients of the time dependent scat- 
tering matrij* 7 -^ are difficult to write in the kink and 
anti-kink basis. 



C. Calculation of the full counting statistics 

The full counting statistics (FCS) is the generating 
function for the moments of the charge transfered be- 
tween two reservoirs through a constriction. In Ref . [30]| . 
the FCS was calculated for the process of opening and 
closing the QPC, with the result that the distribution 
function for the transmitted charge is purely Gaussian. 
Here we shall extend our methods developed in the pre- 
vious sections to find the FCS in the interacting case, 
and conclude that the non-interacting result is in fact ro- 
bust to the presence of interactions. It is plausible that 



V. ENERGY-MOMENTUM NOISE IN THE 
QUANTUM ISING MODEL 

Although in previous sections we showed that the noise 
generated by a quench can also evolve logarithmically 
with time, we will show now that this is not always 
true. This is in contrast to the static entanglement en- 
tropy which depends logarithmically on the size of the 
observed regio n 20 i 60 ' 61 and the dynamic entanglement en- 
tropy which evolves logarithmically with time 42 

Consider a one-dimensional quantum Ising model that 
can be split in two by changing with time the strength of 
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just one link. This effectively implies a change in bound- 
ary conditions at that link. We use the Majorana fermion 
description of the critical quantum Ising model, 



C = ifm{dt - d x )fj R + ir\h{dt - d x )rj L 



(5.1) 



where r\h (vr) is a left-going (right-going) Majorana 
fermion, and f] R (x) = rj R (—x) was flipped. At the center 
of the Ising chain (x = 0) we modify the link strength, 
which translates in the effective Fermion description to 
a coupling to the local energy density operator e = r/LrjR 
at the same point 



£ = i\6(x)r) L r) R 



(5.2) 



This marginal perturbation leads to the following relation 
between the fields on the two sides of the impurity 



with 



M = 




(5.3) 



(5.4) 



When A is fine tuned to the point A = 2, M becomes 
completely off-diagonal M — » ~ia y . When A — > 0, M 
approaches the identity matrix, M — > I. 

The quantum Ising model and its effective field theory, 
the free Majorana fermion, has a (discrete) global Z2 
symmetry. Contrary to the Luttinger-type models which 
have a globally conserved U(l) charge and a locally con- 
served (dimension one) charge current, in the case of the 
quantum Ising model the only global conservation law is 
the energy-momentum and its locally conserved (dimen- 
sion two) energy- momentum current. In contrast with 
the Luttinger model where the quantum quench does not 
change the total charge, in the Ising model the quantum 
quench changes the Hamiltonian and hence the total en- 
ergy. Although the energy-momentum tensor of the Ma- 
jorana fermion field theory (and of the Ising model) re- 
mains locally conserved everywhere except at the location 
of the quantum impurity, these physical differences lead 
to a distinct asymptotic time dependence of the noise in 
the energy-momentum current shown below. 

The energy density along the edge is defined by 



Pe{x) = r\ L id x T] L - fj R id x fjR 



(5.5) 



The rate of back-scattering of energy by the modified link 
is 



Ie = Pb(0 + ) - Pb(P~) 



(5.6) 



Using the relations above, for the case that the link is 
fully transmitting (A = 0), the heat current is zero as the 
contributions on the two sides of the link cancel. In the 
limit of a fully reflecting link, we get pe(0 + ) = — pe{0~), 
and the two contributions add up. 



Since the scaling dimension of the energy current is 2 
we get the following expression for the thermal noise E2'. 



E ? = 7^ let* *i d< 2 ( tl -/ 2 +tj ) 
1 3S 2 At 2 +At 4 
~ ^ 36 2 (6 2 + At 2 ) 2 



(5.7) 



For large At, the thermal noise approaches a non- 
universal, cutoff dependent, value 



1 



1 1 



3ir 2 6 2 tt 2 At 2 



(5.8) 



VI. CONCLUSIONS 

In the previous sections we identified an extended the- 
ory of the bosonic Luttinger description of the quan- 
tum Hall edge states in the presence of an impurity. 
The extended theory, being an j4 m /Za orbifold theory 
(for v = 1/m), explicitly contains an operator inducing 
the transition between the zero-backscattering and fully- 
backscattering fixed points of the impurity. In particular, 
for v = 1 /2 the extended theory decouples into two copies 
of Z2 parafermions, whereas for v = 1/3 it coincides with 
Z4 parafermions instead. The presence of the primary 
field in the algebra which induces the transition between 
the two fixed points at the impurity demonstrates that 
in the sudden approximation equilibrium is maintained 
through each process of opening or closing of the QPC, 
and that no transient effects are expected. 

Using these methods we calculated the FCS for the 
interacting Luttinger liquid for the process of a sudden 
opening (at to) then closing (at t\) of the QPC. Only 
the second cumulant of the transferred charge is nonzero, 
and its logarithmic dependence on t\ — 1 suggests that a 
relation between the entanglement entropy and the FCS 
should exist even in the interacting case. In that sense, 
the noise "measures" the entanglement entropy. The only 
difference from the non-interacting case is the explicit 
appearance of the filling factor in the noise. 

In a more general setting, however, the relation be- 
tween the FCS and the entanglement entropy does not 
seems to be structural. First, it requires the presence of 
a conserved current with dimension 1 (essential for a log- 
arithmic behavior of the correlation function), and of an 
associated and strictly conserved global charge. In con- 
trast, the energy-momentum noise through a weak link 
in the one-dimensional Ising model does not reproduce, 
or even has the same form, as the entanglement entropy. 
We discuss this case briefly in Section [V] Also, for any 
quantum Hall state which contains neutral edge modes, 
the generated entanglement entropy is proportional to 
the total central charge, c, but in the absence of an im- 
purity the noise involves only the pure current present 
in the charge carrying c = 1 theory alone. This can 
be amended by adding the central charge explicitly into 
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the relation between the noise and the entanglement en- 
tropy, but this argument seems to be ignoring the mech- 
anisms that create entanglement in the neutral theory. 
Moreover, even for the simple Luttinger case the central 
charge does not appear in the noise (which depends only 
on the fusion rules and the conformal dimensions of the 
CFT) while it certainly appears explicitly in the entan- 
glement entropy. Indeed, our results show that the noise 
generated is a measurement of dynamical entanglement 
and not of the static entanglement that is associated with 
the observation of a subsystem in an eigenstate of the full 
system42 

Interestingly, we also find that the time dependent cor- 
relation functions have a universal scaling behavior as in 
the case of the global quench.— However, the univer- 
sal behavior is characteristic of an orbifold theory of the 
original system. It may be of interest to study this exper- 
imentally where corrections to scaling are sure to occur. 
The flow away from criticality would presumably involve 
the more exotic fusion rules. 

In summary, in this paper we have computed the noise 
of the tunneling current generated by a quantum quench 
of a point contact in a fractional quantum Hall fluid. This 
was done by explicit solution for the bosonic Laughlin 
state at v = 1/2 and by conformal field theory methods 
for a general Laughlin state. We have also computed the 
growth of the entanglement entropy due to the quantum 
quench. We found that even though the time dependence 
of the noise has the same form as the entanglement en- 
tropy the latter has a more intricate dependence on the 
properties of the conformal field theory and the Hilbert 
space of the edge state. This result suggests that the 
Klich-Levitov protocol may not generally supply a proce- 
dure to measure the entanglement entropy. The question 
of measuring the entanglement entropy remains open. 
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APPENDIX A: DETAILS OF v = 1/2 



so that {i/j(t,x),il)i(t,x')} = S(x - x') holds. The full 
Hamiltonian can now be written as 



dxip {x)id x ip{x) 

2^8(x) [A^(x) 7 + A* 7 ^ t (a;)] (A3) 



Note that p (x) = ^ (x)tp(x), and {7,7} = 2. Here 7 is 
a Klein factor which appears in a careful handling of the 
zero modes of the boson field. 

Assuming for simplicity A = A* G R, the equations of 
motion can be found and are given by 



(A4) 



idtip = id x ip — v2tt\'-/5(x) 
idtip^ = id x ip' + V^nXjS^x) 
id t 7 = 2v / 2^A(?A t (0)-i/>(0)) 



The equations of motion can be solved by expanding the 
fields in modes as (see e.g. I62I ) 



x < 
B^ x > 



til ^ 

= E 



aL„ x<0 
x>0 



(A5) 



The modes satisfy the following relations, derived by in- 
tegrating the equations of motions around x = 0, 



i(Boj — AJ) = V2^A 7w 
i{Bl u - Alj = - V2^A 7a 

uj 1u} = 2V2^A [|(A, + BJ) - \{Al u 
Defining cjb = AttX 2 , we get 



(A6) 



Bl 



Bl 



1 



= M, 



IUJ — UJb 

A t 
Al 



IUJ —UJB 

-ujb iuj 



A, 
Al 



(A7) 



So that when ujb — * (fully transmitting QPC), 
M u — > /, while when ujb — * 00 (fully reflecting QPC) 



In this appendix we give some details of the re- 
fermionization technique used in the exact solution of 
the v — 1/2 (bosonic) quantum Hall effect. For v = 1/2, 
the Lagrangian for the odd boson in the presence of a 
QPC at x — is given by (v = 1) 

C = ^-d x Mdt - d x )<t> + \5{xy*°^ + h.c. (Al) 

47T 

where the odd density p a (t,x) = [pB.{t, x) — 
PL(t, -x)]/v2. We refermionize by defining^ 

ip(t,x) = :e^^ : (A2) 



APPENDIX B: MODULAR 5-MATRIX AND 
FUSION RULES AT LEVEL l/v = m 

With the convention that the matrix takes the vector 
labelled by (1, j, <f>m , <f>k' , o~j , Tj) to the same vector with 
j — > i, k — > k' and cr,j = 25y — 1, the modular S'-matrix 
for m even can be written as: 
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/ 1 
1 

1 

2 



2 2 2(-l) fc 4cos7r 2n , 
\fm —\pm (Jijy/m 
V \M -\fm (Jijy/m 

while for to odd, the modular S'-matrix takes a slightly different form. 



1 1 2 

1 1 2 

1 1 2(-l) fc ' 

2 Of 1 \k A _ fefe' 



Ira 



in 



— (5jj\/2m (5^- 



'to 
-■y/m 





f 2m J 



S. 



rn.odd 



/8m 



1 


1 


l 


2 


\/m 


V™ \ 


1 


1 


l 


2 


— \Jm 


-\frh 


1 


1 


-l 


2(-l) fe ' 


ia t j ^frh 




2 


2 


2(-l) fc 


4cos27r^ 








/to 


-\fm 


iaij y/m 







-e iw(7i ^ 4 \/2m 



(Bl) 



(B2) 



For to even, the elements 4>m} form a Z 2 x Z 2 
subalgebra. One has 



1 



. x ^ = 1 

The twist operators have the fusion rules 

a l x cr t = 1 + + ^ 0fe 
i 

£71 X £72 = ^ 0fc 



(B3) 



k.even 



(B4) 



For to odd, the elements {1, j, </>m } form a Z 4 subalgebra. 



i x 3 = 1 

x = i 



(B5) 



and for the twist fields, 



<ji x cr,; 



(J\ X £7 2 = 1 + 



k.odd 

E * 

k.even 



(B6) 



The fusion rules for the vertex operators <j)k are the same 
for to even or odd. 



4>k x 4> k > = 4>k+k + 4>k-k> 
4>k x 4>k = i + j + 4>2k 

<p m -k X (f> k = fok + <t>m + <t>m 

j X (j) k = 4> k 



(B7) 



The fusion rules for the r» fields can be found easily by 
applying j x tij = t;. 

APPENDIX C: SCHWINGER-KELDYSH AND 
ENTANGLEMENT ENTROPY 



The non-equilibrium problem of locally quenching the 
system at the point contact is best dealt with in terms of 
the density matrix formalism. Initially, the system has 
a density matrix po = IV'oC^)) (V'o (^) I where the states 
\ipo(x)} are eigenstates of the Hamiltonian with Neumann 
boundary conditions imposed at the point contact and 
ends of the system. At time t = 0, the point contact is 
opened and the boundary condition at the boundary is 
changed to Dirichlet. 

After the quench, the density matrix is given by 



J 



W(x)\p(t)\<f>(x)) = {<t>'{x)\ uHt)poU(t) \<f,{x)) = (cy(z)h 

I 



-itH-eH 



\i/>o(x)){il>o(x)\e 



itH-eH 



\</>{x)) 



(CI) 



where the regulator e has been included to adiabatically term describes forward evolution from the initial state to 
cut-off the high momentum modes. Here, H is the Hamil- the final state \4>(x)) while the second describes the time 
tonian with transmitting boundary condition. The first reversed evolution back to the initial state. 
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Following Refs. [40.41], the Schwinger-Keldysh time 
contour can be thought of as the time contour for a con- 
formal field theory. On the forward branch, one continues 
t = — 6 — it, 4>{x, t) — > 4>{x, t). Now, r — * —it, Wick ro- 
tating, this can be thought of as propagating the initial 
state for time t to the final state. 



{<f>'fa)\e 



-UH—eH 



= W(x)\e- lHT \Mx)) (C2) 



One can think of the density matrix in terms of a path 
integral, but to do so, a complete set of states needs 
to be introduced. Here we work with the extended 
Hilbert space which is spanned by states with all possible 
boundary conditions. Inserting a complete sets of states 
/ d<j>i \4>i){4>i\ at each time slice Sri and also a complete 
set of conjugate states J djli iTVi/fa]. 



{cf>{x)\e- lHT \Mx)) 



N N 



»=1 z=l 

Dcj> ef dTdx £ I^' r )l 6(<P(x,t = 0) - 4>{x)) Stefan) - ipofa) 



(C3) 



Now, dr = —idt so that this is the Euclidean path inte- 
gral with the boundary condition that at n = — e — it, 
the field takes value 4>{x, n) = ipofa) and 4>{x, r) = 4>fa) 
at t = 0. A similar set of manipulations can be per- 
formed for the backward branch as well. At the point 
t = 0, they differ on a set of measure zero so that the 
path integral can be thought of as over the entire surface 
parameterized by (x, r) with changing boundary condi- 
tions at n and T2- At the end of the calculation, one 
should continue back to real time t and take the limit 
e — » 042 Here, before the quench t < 0, we apply reflect- 
ing boundary conditions and after the quench, we have 
transmitting boundary conditions. In the sequel, we de- 
note this surface, with changing boundary conditions in 
time at x — as E. 

To compute the entanglement entropy, one uses the 
replica trick to compute, Tr p\. The reduced density ma- 
trix can be visualized as stitching together the cylinder on 
the B side while leaving the A side open. Multiple copies 
of pa are stitched cyclically on the cut i 40 ' 41 In Ref.[2Cj], 
it was shown using general transformation properties of 
the stress energy tensor that this branch cut amounts 
to a twist operator <f> n (z) insertion on the Riemann sur- 
face with conformal dimension, A„ = ^ (n — 1/n) . The 
stitched surface amounts to computing the correlation 
function of the twist operator. 



Tr p\ = ($„(z = 0)> 3 



(C4) 



Here 



2 denotes an average taken over the surface with 
changing boundary conditions. To compute this corre- 



lation function, we use the Zhukowski mapping which 
avoids the algebraic structure discussed in the text. The 
advantage is that including an n-twist field as part of an 
extended Hilbert space is difficult and is left as an open 
question here. Using the mapping, 



W 



-(z+ V* 2 + e 2 ) 



(C5) 



the surface E is mapped to the right half plane with a 
single boundary condition, E'. The correlation function 
on E is related to the correlation function on E' 



dw 

dz 



-A„ 



(C6) 



Differentiating with respect to n at n = 1 yields the en- 
tanglement entropy 



s(t) = £ los (!!±f! 

6 \ ae/2 



ci 



(C7) 



where c\ is a non-universal constant and a a lattice 
cutoff In the limit t ^S> e the entanglement entropy 
becomes 



S(t) = I log i + 

6 a 



(C8) 



so that the amount of entanglement grows logarithmi- 
cally in time after the quench at t = and is proportional 
to the central charge. 
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